LINEAR STATE SPACE THEORY IN THE WHITE 
NOISE SPACE SETTING 
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Abstract. We study state space equations within the white noise 
space setting. A commutative ring of power series in a countable 
number of variables plays an important role. Transfer functions are 
rational functions with coefficients in this commutative ring, and 
are characterized in a number of ways. A major feature in our ap- 
proach is the observation that key characteristics of a linear, time 
invariant, stochastic system are determined by the corresponding 
characteristics associated with the deterministic part of the sys- 
tem, namely its average behavior. 
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In a preceding paper, see [J] , the first two authors began a study of lin- 
ear stochastic systems within the framework of the white noise space. 
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There, the emphasis was on stability theorems associated with convolu- 
tion systems. The present paper is concerned with state space theory. 
Specifically, we study systems defined by state space equations when 
randomness is allowed in the matrices defining these equations. 

Remark 1.1. While the discussion to follow is restricted to discrete 
time, most results apply to continuous time in an obvious way. Con- 
tinuous time is not explicitly pursued in this work. 

To set the problem in perspective and provide motivation, we first recall 
some well known facts from linear system theory. There, state space 
equations of the form 

x n +\ = Ax n + Bu n , n = 0, 1, . . . 
Dn = Cx n + Du n , 

play an important role. In (|TIjl . A G C NxN , B G C Nxq , C G C pxN , 
D G C pxq , the states x n take values in C N , the inputs u n in C q and 
the outputs y n in C p . Taking the Z transform, assuming x = 0, ( 11.11) 
leads to 

j(0 = W() + (W(0 

^ ' ' &(C) = C&(() + D<&(C), 

where the Z transform variable is denoted by C, so that 

(1.3) y(c) = jr(cwc), 

where 

oo oo oo 

n=0 n=0 ra=0 

and 

(1.4) Jf(() =D + (C(I N - CAy'B. 

The function J^(C) is called the transfer function of the system defined 
by (11. ip . It is also possible to replace C by a field & over the complex 
numbers. Matrices A, B, C and D then have their entries in ^, and 
J^iC) makes sense for all ( such that (In — (A) is invertible. The case 
where the components of A, B, C and D belong to some commutative 
ring is of special interest. See for instance [15], [T7j, [16], [11], [12], [5]. 
This ring is often assumed Noetherian, to enable to formulate results. 
Even when the ring contains the complex numbers, formula (11.41) does 
not make sense in general because elements in the ring need not be 
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invertible in the ring. It will make sense in a normed ring for ( small 
enough, as is seen using the power expansion 

oo 

(i.5) (i N - ca)- 1 = c An - 

n=0 

The ring D\ defined below is not a normed ring, but still it will be pos- 
sible to define expansions of the form (jl.5p in it. 

As we have explained in our previous paper pQ, a Gaussian input into 
a linear system with nonrandom coefficients, will result in a Gaussian 
output. In the present series of papers, and in particular in the present 
work, we aim to model linear Gaussian input-output relations when 
the underlying linear system is random. Here we allow Gaussian in- 
puts and randomness in the matrices A, B, C and D, in such a way 
that the outputs remain Gaussian. While indeed a Gaussian input into 
a linear system with random coefficients cannot be expected to result 
in a Gaussian output, we will use the white noise space setting (see [9], 
[TO] and Section [2] below), and replace the pointwise product with the 
Wick product, enabling Gaussian input-output relations when the un- 
derlying system has random coefficients. This framework will preserve 
Gaussian input-output relation, while allowing uncertainty in the form 
of randomness in the linear system under study. Such a setting may 
prove useful so as to study a linear state space system with nonrandom 
uncertainties, a system that indeed maintains Gaussian input-output 
relation. This, by utilizing the Bayesian embedding approach to solve 
problems associated with a system subjected to a deterministic uncer- 
tainty, through solutions to corresponding problems associated with 
systems with random uncertainties, see e.g. [13] . 

In the white noise space setting, the space of complex numbers C is 
replaced by a space of stochastic distributions called the Kondratiev 
space, denoted by S-i. This space contains C, is the inductive limit of 
a certain family of Hilbert spaces (see (12.21) below), and is nuclear; see 
[TO] Definition 2.3.2 (b) p. 30, Lemma 2.8.2 p. 74]. A key element in 
our formulation is a product defined on S-i, namely the Wick product, 
denoted by u()v, which reduces to multiplication by a constant when 
one of the elements u or v is non random. We thus replace the equations 
CHI by 

x n+ i = A0x n + B0u n 

(1 ' 6 ^ Vn = C0x n + D0u n 
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where A e (S^) NxN , B e (5'_i) JVx «, C E (S^) pxN , and D £ (S^) pxq . 
The states x n take values in (S_i) N , the input u n in (S_i) q the output 

in (s^y. 

A fundamental tool in white noise analysis is the Hermite transform 

F i-> 1(F) 

(see below), which associates to every element in SLi, a power series 
in a countable number of complex variables, 

(1.7) z = (zx, z 2 , za, . . .) 

and transforms the Wick product into a point- wise product: 

(1.8) l(FoG)(z) = (l(F)(z))(l(G)(z)), VF, G e 

The image of the Kontradiev space under the Hermite transform is a 
commutative ring without divisors of zeros (that is, a domain), which 
we will denote by £H. It is not Noetherian, so most results in system 
theory on commutative rings cannot be applied. Still, it has a very 
important property, which allows us to proceed. An F e !SR is invert- 
ible in D\ if and only if its constant coefficient is non zero (recall that 
F is a power series). More generally, for p G N, an F 6 y{P x P will 
be invertible in 9^ pxp if and only if the matrix F(0) (which belongs to 
C pxp ) is invertible; see Theorem 13.11 below. This theorem follows from 
a non trivial result on the characterization of the range of the Hermite 
transform, given in [10, Theorem 2.6.11, p. 62]. 

We now take the Z transform and the Hermite transform of (11 .6p . The 
Z transform of the series (I(w n ))n=o,i,... is denoted by ^(C, z) 

oo 

^(c,^) = EniK))w, 

n=0 

and similarly for < 3f((,z) and ^~((,z). We obtain: 

(I N -(A(z))%-(t,z) = B(z)W((,z), 
{ ' } nCz) = C(z)%-((,z)+-D(z)W((,z), 

where 

A(z) = I(A)(z), B(z) = I(B)(z), C(z) = I(C)(z), D(z) = I(D)(z). 

Since A(z) is bounded in a certain neighborhood of the origin (see 
Theorem 13. II below) . the matrix (ijy — (A(z)) is invertible for ( small 
enough, and we can write 
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where 



(1.10) 



Jf(C, z) = D(z) + {C(z){I N - CA(z))- 1 B(z), 



is the transfer function of the system defined by the equations (II. 61) . 

Remark 1.2. When we set z = in (11.91) we retrieve (1 1.4ft . i/iai is, 
we are back within the deterministic setting. 

We now view (jl.lOp as a 9^-valued function. Consider ( such that 



It follows from Theorem 13.31 below that (Jjy — (A) is invertible in d\ 
for such £. Therefore the 9^-valued function given by 



is well defined for £ satisfying (II. lip . 

Functions of the form fll. 12j) will be called rational functions associated 
with the white noise space. We note that in [2] another approach to 
rational functions, with emphasis on rationality with respect to a finite 
number of the variables Zk is considered. The purpose of this paper 
is to give a number of equivalent characterizations of functions of the 
form (11.121) and to study the notions of controllability, observability 
and minimality in the setting of the ring !>K. 

A major non-trivial feature of this work is the observation that key 
characteristics of a linear, time invariant, stochastic system (e.g. in- 
vertibility), are determined by the corresponding characteristics associ- 
ated with the deterministic part of the system under study, namely its 
average behavior. For instance, a realization of the perturbed system 
(I1.12p will be observable (see Definition 15 . II below) if the corresponding 
realization of the unperturbed system, namely with z = is observable. 
See Theorem 15.31 

The paper consists of six sections besides the introduction, and its out- 
line is as follows. We review, as already mentioned, white noise space 
theory in Section [2j We study the ring 9^ in Section [31 Equivalent char- 
acterizations of rational functions are given in Section HI Observable 
pairs are studied in Section [6j In Section 6 we consider controllable 



(1.11) 



det(I N - CA(0)) + 0. 




^(C)=D + CC(/ iV -CA)- 1 B 
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pairs, and briefly discuss minimal realizations. The last section con- 
siders the case where the functions take values in one of the Hilbert 
spaces which make 9^. 



2. A BRIEF SURVEY OF WHITE NOISE SPACE ANALYSIS 

The starting point to construct the white noise space is the Schwartz 
space S of real-valued smooth functions which, together with their 
derivatives, decrease rapidly to zero at infinity. For s G S, let 
denote its L^M) norm. The function 



is positive (in the sense of reproducing kernels) for s±, s 2 running in S. 
The space is nuclear. By an extension of Bochner's theorem to nuclear 
spaces due to Minlos (see [H], [H Theoreme 3, p. 311]), there exists a 
probability measure P on S' such that 



where we have denoted by (s 1 , s) the duality between S and S'. The 
real Hilbert space L 2 («S', J 7 , dP), where T is the Borelian a-algebra, is 
called the white noise space. We will denote it by W, and its elements 
by u, by setting Q = S f . 

Among all orthogonal Hilbert bases of the white noise space, one plays 
a special role. It is constructed in terms of Hermite functions and its 
elements are denoted by H a , where the index a runs through the set i 
of sequences (a±, «2, • • •), whose entries are in 



and ttfc 7^ for only but a finite number of indices k. Furthermore, 
with the multi-index notation 



K{s x - s 2 ) = e 




N = {0,1,2,3,...}, 



we have 




H, 



|2 



a\. 



a\\W ~ 
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extends to a unitary map between W and the reproducing kernel 
Hilbert space with reproducing kernel 

k(z,w) = e {z ' w)e i = 
where z, w run through £ 2 - 



The Wick product in W is defined by the formula 
Ha<)Hp = H a+ p, a, ft G t, 
and the Hermite transform is defined through linearity as 

(I(ff a ))(s)=s". 

The space W is too small to be stable under the Wick product, and 
one defines the Kondratiev space S-i, within which the Wick product 
is stable. More precisely, S_i is a nuclear space, and is defined as 
the inductive limit of the increasing family of Hilbert spaces Tik, k = 
1,2,... of formal series Ylaee faH a such that 

(2.2) \\fh d = (j2\f«\ 2 ( m r ka ) <°°' 



where, for ft G £, 

(2Nf = 2 ft (2x2) fc (2x3) ft -- - . 

That the Wick product is stable within S-i is made more precise by 
Vage's inequality (see [101 Proposition 3.3.2, p. 118]), which we now 
recall. Let / and k be natural numbers such that k > l+l. Let h e Tii 
and u G TCk- Then, 

(2.3) \\hOu\\ k < A(A;-Z)||/i||,|H| fc , 
where 

(2.4) A(k-l) = J2( m ) {l ~ k)a - 

For a proof that A(k — I) is finite, see [TQl Proposition 2.3.3, p. 31]. 
The series J2 a ee s <± za ^ e sa ^ ^° ^ e convergent at z if 



E 



z < oo, 



that is, if it is absolutely convergent, see [10, p. 60]. The following 
easy lemma will be used below. 
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Lemma 2.1. Assume that f(z) = Y. a& tfaZ a and g(z) = J2 a ee9aZ a 
are absolutely convergent power series at z. Then 

(2.5) i/(*M*)i<5>r-l Yl /^I<(Ei^ii z H(Ei»«in°)' 



76i? a+/3="f 



and in particular, the product fg is an absolutely convergent power 
series at z, and it holds that 

( \ 



f(z)g(z) = 



7&? 



E /« 



90 



«+/3=7 
\ ot,Pe£ / 



Furthermore, (f(z)) n is an absolutely convergent power series at z for 
all fiGN. 



Proof: The power series 



5> 

7Q? 

is absolutely convergent since 



/ 



\ 



«+/3=7 



/ 



E 

7e£ 



q+/3=7 



7^ 



o+/3=7 



(£i/«ii*n(£i0«iw a )- 



aee 



□ 



3. The ring I(5_i) 

Consider the image 9^ d =' I(iS'-i) under the Hermite transform of the 
Kondratiev space. This is a space of power series which has been 
characterized in [TQl Theorem 2.6.11, p. 62]. In that statement, (C N ) C 
denotes the space of finite sequences of complex numbers indexed by 
the integers, and the set K q {5) is defined by 

(3.1) K q (5) = {z e C N : k1 2 (2N) 9Q < 6 2 }. 

Q^(0,0,-) 
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Note that a = (0, 0, . . .) is excluded from the sum. See pm Definition 
2.6.4, p. 59]. 

Theorem 3.1. [10, Theorem 2.6.11, p. 62] 

(1) If F(u) = ^a^M G S-i, then there exist q < oo, M q < 
oo such that 

(3.2) < K\ H < M q (j2( m ) qa \ z ' 

for all z G (C N ) C . In particular, 1(F) is a bounded analytic 
function on K g (5) for all 5 < oo. 

(2) Conversely, suppose g(z) = b a z a is a given power series of 
z G (C N ) C with b a G C, with q < oo and 5 > such that g(z) is 
absolutely convergent for z G K q (5) and 

sup \g(z)\ < oo. 

z€K q (8) 

Then there exists a unique G G S-i such that 1(G) = g, namely 
G(u) = J2 b *H a (uj). 

A characterization of convergent sequences in S-i is given in the fol- 
lowing theorem proved in [10] . It will be used in particular in the proof 
of Proposition 15.21 

Theorem 3.2. [10, Theorem 2.8.1, p. 74] A sequence of elements F^ 
in the Kondratiev space S^i converges to F G 5_x if there exist 5 > 
and q < oo such that l(F^ n ') converges to 1(F) pointwise boundedly, or 
equivalently, uniformly, inK q (8). 

The main result of this section is: 

Theorem 3.3. d\ is a commutative ring, which contains C and has no 
divisors of zero. Furthermore, let x(t) = Y^=o Xn ^ n be a P ower series, 
with strictly positive radius of convergence. Let p G N. Then for every 
r G W xp such that r(0) = pxp , the series 

oo 

(s(r))(z) = 5> B (r(z))» 

n=0 

converges to a limit in 9F xp . Ify(t) = YlnLo ^ s an °ther such power 
series, then 

(3.3) (xy)(r) = x(r)y(r), Vr G 9t. 
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In particular, an element s is invertible in 9l pxp if and only if s(0) is 
invertible. 

Proof: To simplify the notation we give a proof for p — 1. The fact 
that we have a ring follows from the formula (jl.8p . The way to prove 
the second claim is to use Theorem 13. II to show that the a-priori formal 
power series 

oo 

E*»(r(*)) B 

n=0 

is in fact the image under the Hermite transform of an element in S-±. 
Since r is the image of an element of 5_i under the Hermite transform, 
it satisfies ( 13. 2\\ for some q G N and a constant M q > 0. Let r x be the 
radius of convergence of the power series defining x (and similarly for 
r y below). We choose 5 such that 

M q 5 d = p <r x , 
Then, by (13.2j) . we have for z G K q (5), 

H z )\ < p, 

hence 

oo oo 

\J2 x n(r(z)) n \<J2\ x n\p n , zeK q {5). 

n=l n=l 

We conclude the proof by using Theorem 13. II We first prove (13.31) . By 
the preceding arguments we know that x(r), y(r), and (xy)(r), are well 
defined. On the other hand, for \r(z)\ < min (r x ,r y ) we have: 

oo n 

(x(r(z))(y(r(z)) = x p y n ^)(r(z)) n 

n=Q p=Q 

We now turn to the last statement. Assume that s is invertible in ISH, 
and let u G R be such that su = 1. Then, in particular, s(0)u(0) = 1, 
so s(0) 7^ 0. Conversely, we can assume without loss of generality, that 
s(0) = 1. It suffices then to take in (Q x(t) =l-t, y{t) = (1 - 1)" 1 , 
and r = u - 1. (Note that r(0) = 0.) □ 

4. Rational functions 

Let f(C,z) = Yl^Lo fn( z )( n ^ !LH pxp (C) be a power series with coeffi- 
cients in DV xq . Define 

*o/(c.*) = /(Cz) - /(M 
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Theorem 4.1. LetH{() = Er=o f «C e W x<? ((0)) be a formal power 
series. Then the following are equivalent: 

(1) Components ofH are obtained by adding, multiplying and dividing 
polynomials of ^[Q, with division being performed only when the con- 
stant coefficient is invertible in 91. 

(2) Tt admits a realization in the form of (11.121) . with coefficients ma- 
trices having entries in 9^. 

(3) The formal power series converges in a neighborhood of the origin, 
and there exists a finite number M such that for every n > M the func- 
tion RqH is a linear combination of R H, . . . Rq I ~ 1 H with coefficients 
in ">R. 

Proof: We first note that elements of the form (11.121) are convergent 
power series in R and not only formal power series, as follows from 
Theorem 13.31 Elements 7i G 9\ pxq (() of the form 

(4.4) W(C) = D or W(C)=CC, 

where C,D G 9F X|J are clearly in the form (I1.12p . Furthermore, as is 
well known, if p = q and 

= D + CC(Jn - CA) _1 B G W xp (C) 

with D invertible, then we have: 

(Jf (C)r 1 = D 1 - CD- 1 C(J JV - CA X )~ 1 BD- 1 , 

where 

A x = A-BD _1 C. 

Furthermore, if 

M[{C) = Di + CCxCJ^ - CAi)-^! G W^ S (C) 

and 

jr 2 (c) = d 2 + cc 2 (/,v 2 - ca 2 )~ 1 b 2 g ^ SX92 (C), 

then 

(^!^)(C) = D + CC(/tv - CA) _1 B, N = N, + N 2 , 
with D = D!D 2 and 

A -( A o' B i?)- B -( B ^ 2 )' c ^ (Cl DlC2) ' 

A sum of matrices is a special case of a product, as follows from the 
formula 



M 1 + M 2 =( Ml I p )(£) 



where Mi and M 2 G W xq . 
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See for instance [3j for some of these formulas when the coefficients are 
complex. It follows from these formulas that any 9^-valued function 
(that is, when p = q = 1) which is obtained by addition, multiplication 
and, when defined, inversion, of functions of the form (14.41) . is of the 
form (I1.12p . The matrix- valued case is obtained by concatenation using 
the formulas 

(Mi je 2 )(o = 

= (D, D 2 ) +C ( Cl C 2 )( Wa -c(^ a,))"^ 1 B 2 )> 

and 

(5) «) - 

= (d 2 ) +c ( < o 1 cj - c {o a 2 )) (b 2 ) 

for two functions M{ and of appropriate dimensions which admit a 
realization. Thus, (1) implies (2). We now prove that (2) implies (3). 
Let ^ be of the form (JTI2"]) . Then, 

i^Jf (C) = C(Jjv - CA^A^B, n = 1, 2, . . . 

But the Cayley-Hamilton theorem holds in any commutative ring (see 
for instance p§ p. 14], [6, Theorem 4.3, p. 120], | p. A III.107]). 
Therefore there exist an M e N, a monic polynomial p of degree M 
with coefficients in such that p(A) = 0. It follows that for n > M, 
the function R^Ti is a linear combination of 1, RqH, . . . R^^Ti with 
coefficients in !5R. 

We now assume that (3) is in force and prove that (1) holds. First, 
assume that TL is 91- valued (as opposed to l >R pxq - valued). By hypothesis 
there exists a matrix A G y{ MxM such that 

R (i R n ■ ■ ■ R^- X n) = (i r h ■ ■ ■ rm- x h) a. 

Hence 

(1 RqH ■■■ RM- l H) = (l RoH ■■■ R^-^miM-CA)- 1 . 
Thus 

R H = C(I M -(A)- 1 B, 
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with 



C = (1 R H ■ ■ ■ R^H) (0) and B 



1 



w 



and the result follows. The matrix-valued case is treated in much the 
same way. □ 



5. Observable pairs 

Consider an 9^-valued function of the complex variable (, of the form 
dTTJ): 

jr(C) = D + CC(/ w -CA)- 1 B. 

By Theorem I3.3| we know that Jf? is well defined, in particular for £ 
such that 

det(I N - CA(0)) ? 0. 

Setting z = in (I1.12p we get the unperturbed transfer function, which 
motivates Theorem 15.31 below. We first give a definition and a propo- 
sition. 

Definition 5.1. The pair (C, A) G W xN x 9K NxN is called observable 
if the map 

f i-> (Cf CAf CA 2 f • • •) 

is injective from £H n into (W) N . 
See [3 §2.2 p. 58]. 

Equivalent ly we have: 

Proposition 5.2. Realization (II. 12ft is observable if and only if (with 

fe^j 

C(I N - CA)- X f = 0j£" f = o£. 

Proof: By Theorem 13.31 with x(t) = (1 — C^) -1 j we have: 

oo 

(I N -(A)- 1 = Y,( P & P - 

p=0 
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Therefore, and using Theorem 13.21 

oo 

C(/ iV -CA)- 1 = C^C P A^ 

n=0 

oo 
n=0 

□ 

Theorem 5.3. Assume that the realization 

(5.1) JT(C, 0) = D(0) + CC(0)(/jv - CACO))-^^) 
is observable. Then realization (j!.12p zs observable. 

Proof: Assume first that realization (15.1 ft is observable, and let 
be such that 

(5.2) C(I N - CA)- X f = 0. 

To prove that the realization (I1.12p is observable we need to show that 
all coefficients f a in the expansion f(z) = ^2 a ^ifaZ a , are identically 
zero. Since ( 15. lft is assumed observable, setting z = in ( 15.2ft leads 
to fo = 0. Let us now put z = (z\, z%, ■ ■ ■ , z p , 0, 0, . . .) in (15.21) with 
p > 1 and differentiate with respect to z\. We obtain (with ' denoting 
differentiation with respect to Z\) 

C(z)(I N - CAC*))- 1 ^) + C(z) ((I N - (A(z)y 1 )' f (z) + 

+C(z)(I N -(A(z))- 1 f'(z)=0. 
Setting z\ = z 2 = ■ ■ ■ = z p = we obtain that 

C(0)(/ jv -CA(0))- 1 /(iao,...) = 0, 

and hence /(i,o,q,...) = since the pair (C(0), A(0)) is observable. Dif- 
ferentiating in turn (15. 3p with respect to z±, we obtain an expression of 
the form 

(5.4) X(z) + C(z)(I N - CAC*))- 1 /"^) = 0, 

where X is a finite sum of the form 

M 

X(*) = £lL>)/(*%) 

3=1 
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where r 3 - G {0, 1} and Uj-(z) is analytic in z\ and may depend on (. 
Setting z\ = in f)5.4p and taking into account that f = /(i,o,o,...) = 0, 
we obtain that 

C(0)(/ jv -CA(0))- 1 /(2ao,..)=0, 

and hence /(2,o,o,-) = 0- More generally, an easy induction argument 
shows that the cui-th derivative of (15.21) is of the form 

(5.5) X(z) + C(z)(I N - CA(z))-V K) (z) = 0, 
where X is of the form 

M 

x(z) = ]Tu^)f("%), 

Uj being analytic in and rij G {0, . . . , a\ — 1}. Setting z\ = in 
( 15. 5p we obtain that /( aii o,o,...) = 0. 

Similarly, by setting z = (0, z 2 , 0, . . .), and more generally 

z = (0,0,0,. ..,^-,0,...) 

in (15. 2p . and differentiating, we obtain that f a — for all a G £ which 
have only one non-zero component. We now prove that /n i,o,o,...) = 0- 
To that end, set z = (^i, 2; 2 , -2 3 , • • • , z p , 0, . . .), with p > 2, in (15.21) and 
differentiate this equation with respect to z\ and z-i. We obtain an 
equation of the form 

(5.6) X(z) + C(z)(I N - CAiz))- 1 -^-^) = 0, 

where now X is a finite sum of elements of the form XJ(z)f(z) and 
XJ(z)q£:(z), with j G {1, 2} and U analytic in z\, z 2 , . . . , z p . The fact 
that 

fo = /(i,o,o,-) = /(o,i,o,...) = 
implies that X(0) = 0. Setting z\ = z 2 = in ( 15. 6ft then leads to 

C(0)(/ Ar -CA(0))- 1 /(i,i,o,..) = 0, 

and hence /(i,i,o,...) = 0, where we have used the observability of the 
pair (C(0), A(0)). By successive differentiation and setting z = we 
obtain that /( ai ,a 2 ,o,o,-) = ^ ^ or ever Y choice of natural integers a\ and 
a 2 - The fact that all other coordinates f a are zero, and hence that the 
pair (C, A) is observable, is shown by induction as follows: 

Induction hypothesis: For N G N, it holds that 

(5.7) /(ai,O2,...,aw,0,0,..0 = 0, V(«i, a 2 , • • • , o>n) G (Nq)^, 
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and 
(5.8) 



f)a\-\ hajv 

C(z)(I N -CA(z))- 1 



dz^dz* 2 ■■■dz aN 



N 

Qa\-\ ha N 

= X aij ..., ajv (z) + C{z){I N - (A(z))- 1 dz a ldz a 2 ; .. dz o<N f ( z )> 

where 

(5.9) z = (z\, Z2, ■ ■ ■ , z p , 0, 0, . . .) with p>N, 
and X ait _ }0 , N (z) is of the form 

(5.10) X ai ,..,aM = fl0 ) fl0 ) ^ 



i=i ftz? 1 <9zf 2 • • • dz 



N 



where pf < oti for i = 1, . . . , N and 

(5.11) 13? + ■ ■ ■ + I3 { ff <a 1 + --- + a JV , 

with the functions XJj analytic in the variables Z\, . . . ,z p . 

The induction hypothesis holds for N = 1, as we have shown above. 
Assume that it holds at rank N. We take p > N + 1 in (15.91) and 
differentiate (15.81) with respect to z N+1 . Since 



U iW I p ,0) ^ f W 

dz/ <9,2 2 2 • ■ ■ dz^ z N+ i 



the term 



is of the form 



d"X-ai,...,a N \ z ) 



dz 



N+l 



i=i dz? 1 dz% 2 ■■■dz P N N dz N+1 
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where P is possibly different from M above, the V, are analytic in the 
variables z\, . . . , z p , and the (3^ are as above. Differentiating the term 

aaiH hctjv 

in (15.81) with respect to z N+ i we obtain a sum of two terms. The first 
is, 

while the second takes the form 

f)aiA h«iv+l 

This proves (15. 8p for (an, . . . , a/v, 1), with X Qlj ... )CW) i being the sum of 
fl5TT3|) and of fl5~T2|) . that is 

(5.14) 

f)aiA h«jv + l 

-C(z)(I N - (A(z))- 1 = 



dz^dz?---dz% N dz N+1 

naiH hctjv+1 

= x ° + C(2)(/ " " CA(z)) " <-^...fe7^ + . f(2) 

Setting zi = • • • = z p = in that expression, we obtain 

C(0)(I N - CA(0))- 1 /(a 1 ,..., aiV ,i,o,o...) = 0, 

and hence f( ai ,...,a N ,i,o,o...) — 0. Differentiating (I5.8P a finite number of 
times with respect to -Z/v+ij a similar argument will show that 

f{ ai ,...,a N , aN+1 fl,0...) = 0, VftAT+i G N. 

□ 



Remark 5.4. X7ie converse to the previous theorem does not hold. 
That is, the observability of the pair (C, A) does not imply the observ- 
ability of the pair (C(0), A(0)). As an example, take N = 1 and 

C(z) = z x , A(z) = 1. 

The pair (C, A) is observable since, for f e 91, 

-^—i(z) = => f = 0. 

l-C 

But the pair 

(C(0),A(0)) = (0,1) 

is not observable. 
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6. Controllable pairs and minimal realizations 

In this section we study controllable pairs and minimal realizations 
within the setting of the ring We first recall that, given a commuta- 
tive ring R, one of the characterization for a pair (A, B) 6 ji NxN x R Nx i 
of matrices to be controllable (or reachable) is that the columns of the 
matrix 



generate R N ; see [HI p. 55]. 

In the classical case (that is, for the complex numbers, or more gen- 
erally, for the case of a field), it is well known that the pair (C,A) 
is observable if and only if the pair (A T , C T ) is controllable (with T 
denoting transpose). This duality principle does not hold in general in 
the case of an arbitrary commutative ring. Only the following direction 
holds: 

Theorem 6.1. [5, Theorem 2.7, p. 59] Let R be a commutative ring 
and let (C, A) e R pxN x R NxN . Assume the pair (A T , C T ) controllable. 
Then the pair (C, A) is observable. 

As explained in [51 p. 59], the lack of duality comes form the fact that 
an homomorphism of modules (say /, from the -R-module Mi into the 
-R-module M-i) can be injective without being residually injective. Re- 
call that residual injectivity means that, for every maximal ideal / of 
R, the induced map from Mi/IR into M 2 /IR is injective when / is 
injective. 

Theorem 16.11 does not help us to study controllability based on observ- 
ability. Furthermore, in [HI Theorem 2.3 p. 178], it is shown that a 
necessary and sufficient condition on a commutative ring for the dual- 
ity principle to hold for all pairs is that every finitely generated faithful 
ideal of the ring contains a unit. As a corollary, the authors of [8] state: 



Proposition 6.2. [8j Corollary 2.4 p. 179] J. If the duality principle 
holds in a commutative ring, then the ring is a total quotient ring. 

For the purpose of the present paper we do not need to recall the defi- 
nition of a faithful ideal (see [HI Theorem 1.5 (ii), p. 177]). The total 
quotient ring of a commutative ring R is the set of formal fractions 
associated with the set of elements of R which are not divisors of zero; 
see [5j p. 35]. Thus, in the case of a ring without divisors of zero (as is 
the case for the ring 9t) the total quotient ring is equal to the quotient 
field associated with the ring; see [5j p. 35]. Since 91 is not a field, it 



(B AB ■■■ A N - 1 B) 
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follows that the duality principle is not satisfied on it. 

After these general preliminaries, let us study controllability and min- 
imality in the setting of the ring 91. Let us repeat the definition of 
controllability: The pair (A, B) £ m NxN x m Nxq is said to be control- 
lable (or reachable) if the columns of the matrix 

(B AB •■■ A^B) 

generate 9t N . See [5l p. 55]. We therefore have: 

Proposition 6.3. Assume the pair (A, B) £ yt N >< N x y\ Nx i to be con- 
trollable. Then the pair (A(0),B(0)) £ C NxN x C Nxq is controllable. 

Proof: Since C 91 , for every / £ there exists a £ H Nq such 
that 

/= (B AB ■•• A^^B) (z)a(z). 

Setting z = in this equality we get the controllability of the pair 
(A(0),B(0)). □ 



Remark 6.4. We note the difference between Theorem \5.^\ and Propo- 
sition \6.3l In the former, observability at z = implies observability 
in 9^. In the latter, controllability in 91 implies controllability at z = 0. 

The converse of Proposition 16.31 would be an analogue of Theorem 15.31 
for the case of controllable pairs. But this is not possible for the ring 91, 
in view of Proposition 16.2^ since 9^1 is different from its total quotient 
ring (which is in fact its quotient field since 9^1 has no divisors of zero). 

Still, we can give a counterpart of Theorem 15.31 for controllable and 
minimal realizations with the following ad-hoc definitions: 

Definition 6.5. Realization (I1.12p will be called 9V controllable if the 
following condition holds: Let f £ 9^ lxAr . Then: 

f(I N - CA)- X B = 0^ q f = 0^ N . 

A realization will be called ^-minimal if it is both observable and 9t- 
controllable. 

We can then state: 

Theorem 6.6. Assume that realization fl5.ll) 

0) = D(0) + (C{0)(I N - CACO))"^) 

is controllable (resp. minimal). Then realization (11.121) is 91- controllable 
(resp. ^Pi-minimal). 
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Proof: The first statement is proved as Theorem 15.31 The second 
statement follows then from the definition of minimality. □ 



7. HlLBERT-SPACE VALUED TRANSFER FUNCTIONS 

We recall that the Hilbert spaces Tik have been defined above by the 
finiteness of the norm (12.21) . By I(7ik) we mean the image of Tik under 
the Hermite transform. We note that 

SK = U- =1 I(?4). 

Theorem 7.1. Let 7i be given by realization (jl.l2p . and let l,k be 

natural numbers such that k > I + 1. Assume that, in the state space 
equations (11.61) . the entries of A and C are in Hi and the entries of B 
and D are in Ti^. Then, the transfer function is \(7ik) -valued. 

Proof: Inequality (12.31) expresses the fact that the multiplication op- 
erator 

Th : u I— ► h()u 

is a bounded map from the Hilbert space Hk into itself. Therefore the 
entries of the y\ pxq - valued function CA n B are in 1(7"^). To conclude 
the proof, it remains to show that for every complex number ( such 
that (In — (A.) is invertible, the power series 

oo 

C"CA n B 

n=0 

converges in I(7ik) to J^(Q- But this is a consequence of Theorem 

□ 

Using once more Vage's inequality (12.31) we have: 

Corollary 7.2. Let now m > k + 1, where k is as in the previous 
theorem. Then, the operator of multiplication by Jrf? sends I(7i m )- 
valued signals into l(7i m ) -valued signals. 
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